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We examine the spatial density within extended objects of the momentum component p^, and 
find relativistically exact connections to Fourier transforms of gravitational form factors. We apply 
these results to obtain semi-empirical momentum density distributions within nucleons, and similar 
distributions for spin-1 objects based on theoretical results from the AdS/QCD correspondence. We 
find that the momentum density in the transverse plane is more compact than the charge density. 
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Mapping the distribution of matter inside the nucleon 
is a goal of hadron structure physics. For example, one 
can relate the density of charged and magnetically po- 
larized material inside a nucleon to Fourier transforms of 
charge and magnetic form factors. This old area has seen 
interesting and notable improvement lately, by way of 
finding relativistically exact relations between form fac- 
tor data and two-dimensional line-of-sight projections of 
the charge density and polarization density [H, d, Q ■ 

In a related fashion, as we study here, one can map out 
the distribution of momentum within a hadron. Here too 
there is old information, in particular that (at some res- 
olution scale) the momentum in a hadron is carried in 
about equal measures by quarks and giuons. More pre- 
cisely, one reports the fraction x of the light-front lon- 



gitudinal momentum, p 



.+ — 



p , that is carried by 



individual constituents. In an extended object, one can 
further ask how the total momentum is distributed point 
by point, either in three-dimensions, or as we shall do 
here, in a two-dimensional line-of-sight projection. 

We will study the spatial distribution of the momen- 
tum component finding relativistically exact connec- 
tions to Fourier transforms of gravitational form factors. 
Gravitational form factors are directly related to empir- 
ical data, in that they are matrix elements of the stress 
(or energy-momentum) tensor which can be obtained 
as second Mellin moments of generalized parton distri- 
butions (GPDs) 0, i]. The GPDs are accessible from 
deeply virtual Compton scattering experiments and can 
be loosely defined as the amplitude for replacing a par- 
ton in a hadron with one of a different light-front mo- 
mentum. In the forward limit, GPDs reduce to ordinary 
parton distribution functions, and the first Mellin mo- 
ments give the electromagnetic form factors; for reviews, 
see e.g. [1, 0, i, i, E El • 

In a light-front formalism, one can kinematically con- 
nect the wave functions of any hadrons of relevant mo- 
mentum. Conventionally, the longitudinal direction is 
the z-direction, chosen to point along the vector direc- 
tion ofp = {pi+p2)/2, where pi{p2) is the momentum of 
the initial(final) hadron. Further, one chooses the frame 
so that the momentum transfer q = P2 — Pi has = 0, 
and satisfies q^ = — (j^ = —Q^. Within the light-front 
formalism, one can project the probability density into 



transverse plane I4I and show that it is obtained from 
the 2-dimensional Fourier transform of a form factor. 

Transverse charge densities have been calculated for 
nucleons P, [3, [1] and dcuterons , using empirical elec- 
tromagnetic form factors. Similar to the charge density, 
one can observe that the r++ component of the stress or 
energy-momentum tensor is the density corresponding to 
the p"*" component of the 4-momentum, 



P+ = / T++ (fxj_dx+ 



(1) 



In the transverse plane in position space, the density 
of momentum p"*" can be defined from expectation value 
of the local operator T~^^{x~^ ,x~ ,x±), 



1 



2p^ 



*|T++(0, 0,6^)1* 



(2) 



such that / (Pb±p'^ {b±) = 1. The state j^*) is localized 
with its transverse center of momentum located at the 
origin, and is formed by linear superposition of light- front 
hehcity eigenstates [p"*", p_L, A), 



, i?^ = 0, A) = 



(2^)2 



P+,P±,X)- (3) 



Normalization factor TV satisfies \Af\'^ J d^p^/ {2it)'^ = 1. 

The matrix elements of the stress tensor between mo- 
mentum eigenstates can be written as 



r++(o) 



xr+,,(Q2), (4) 



where Ai(A2) denotes the initial (final) light-front helic- 



ity, and q± = Q{cos(f>qCx + sinc^ 



The p'^ momentum 



density can be written in terms of 2-dimensional Fourier 
transform of the form factor 



fqi 



(5) 



We will calculate the spatial distribution of p+ in the 
transverse plane for spin- 1/2 and spin-1 hadrons, specif- 
ically the nucleons and the rho mesons. 
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The momentum density results for nucleons will be 
semi-empirical. We use matrix elements of the stress ten- 
sor obtained from second moments of GPDs [1]. "Semi" 
in semi-empirical is a reminder that the GPDs are not 
yet measured in detail. However, the models are con- 
strained to accurately represent measured parton distri- 
bution functions in one limit and to give measured elec- 
tromagnetic form factors in another. 

For spin-1/2 particles, matrix elements of the stress 
tensor involve three gravitational form factors. 



(6) 



2m 



where j^f^p"') = {'yt^p" +p''jt^)/2 and pi,2 =pTq/2. 

For the vector GPDs for spin-1/2 particles, the quarkic 
contribution is defined by 



dy 
2ti 



u{p2,\2) [H'^ix, 1 1) r + E'^ix, t) ^^^"qA uipi,\i)- 

(7) 

(In the frame we use, ^ = — 2q+/p+ = 0.) We also have 
and due to gluons. From the definitions, 

MQ') = J xdxH-{x,^,Q'), 



(8) 



a=g.q 



For the GPDs we use the "modified Regge model" of 
Ref. [ll]. There we find, 



H^{x,0,Q'' 



„q'(1-x)Q= 



qvix)x° 



with 



a' = 1.105 GeV" 
Ku = 1.673, 
Kd = -2.033, 



r]u = 1.713, 
= 1.475, 
Nd = 0.9083. 



rjd = 0.566, 



(10) 



The values of k"''^ are obtained from Kp = (2/3)k„ — 
(l/3)Kd and k„ = -(1/3)k„ -|- (2/3)/vd. 

For the momentum density we need sea and gluon dis- 
tributions. The gluon GPD H3{x,0,Q'^) is discussed 
in [l3|. For applications, they give two forms, each of 
which is a gluon parton distribution multiplied by either 
a dipolc form or a gaussian form for the dependence. 
The results from these two forms are very close to each 



other. A third possibility is to use the same depen- 
dence as for the valence quarks, Eq. This gives a 
result in between the results from the two choices in [31 , 
so will be the choice we use in this paper. We follow a 
similar procedure for the sea quark contributions. 

All quark and gluon distributions are taken from the 
MRST2002 global NNLO fit at input scale = 1 
GcV2, 

u^{x) = 0.262a;""-^9(l - xf-^° (l + 3.83x°-^ -I- 37.65x) , 
d^x) = 0.061x-°-S5(l - x f°^ (1 + 49.05a;°-^ + 8.65a;) , 
S{x) = 0.759x-i-i2(l - xf-^^{l - 1.342;°-^ + 7.40a;), 
g{x) = 0.669x^1(1 - xf-^^{l + 6.98x°-^ - 3.63x) 

-0.23a;-i-27(l-x)«^ (11) 



With these definitions, we calculate 



), (12) 



where, Jq is the Bessel function, b = \bj_\, and Q = |<7j_|. 

Results for p^^^{b) are shown in Fig.[l] both as a con- 
tour plot and along a line through the center of the nu- 
cleon. The result is the same for the proton and neutron, 
assuming isospin symmetry. For comparison, the charge 
density of the proton is also shown. 

The RMS radius of the p~^ density is 0.61 fm, notably 
smaller than the Dirac radius of 0.80 fm (corresponding 
to charge radius 0.87 fm) in this paramctrization. 

Nucleons can also display density shifts due to polar- 
ization. We polarize the nucleons in the transverse plane, 
and denote the polarization direction by S± = cos (^x + 
sin 05 e.y. The state with spin projection = 1/2 in this 
direction will be \p'^ , q±/2, s±} . The impact parameter 
direction is denoted by b± = cos (pb + sin By. The p+ 
density of the transversely polarized state is 



1 f d^q± 



2p+'J (2^) 



(13) 



p+{b^) + sin(0, - 0s) r ^ ^JiibQ) B{Q^). 
1 Jq 2tt 2m 



For the valence quarks, B{0) w —0.02, and p^^^ is es- 
sentially the same as p^^2 1 contrast to the correspond- 
ing electromagnetic case [11,131, and we show no plots. 

For spin-1 particles, the stress tensor matrix elements 
can be expanded in term of six Lorcntz structures multi- 
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FIG. 1: Upper panel: Momentum density (for "p'^") of the 
nucleon, p]y2' projected onto the transverse plane. Lower 
panel: The solid red line is p^^j for the nucleon, along an axis 
through the nucleon's center. For comparison, the proton 
charge density is also shown, as a dashed green line. The 
charge density spreads out more than the momentum density. 



plying six gravitational form factors, 

(14) 



1 



D{q') 



)/3 



where aH^^ = {a'^b'^ -a^b"')/2. 

The two independent helicity conserving form factors 
T+, and To+, defined by Eq. Q, can be expressed in 
terms of the above form factors, 

r+ = A + 7^E^ Z2iQ^), (15) 
T+ = {{1 + 2tj)A + 47^B-2ijD-2tj^E) = ZiiQ^), 

where rj = /{Arn?). The first equality in each line 
above is generic; the second uses an AdS/QCD model 



worked out for rho (and other) mesons in [l^ , where 
1 r dz 



Z2{Q^) = I ^e-^n{Q,z)^{z)^{z) 



dz 



(16) 



We obtain the densities 



^QJ,{bQ)Z2{Q') 
Zn 



dQ 

2lT 



QMbQ)ZiiQ'- 



(17) 



In the hard-wall model [20|, $(z) = and the z- 
integration in Eq. (flBl) is from to zq- The profile func- 
tion Ti.{Q, z) and wave functions ipi^) £^re 



n{Q,z) 



fz^ 



KliQzo) 



V2 



zoJiirrinZo 



h {Qzo 

zJi(mnz) 



l2iQz) + K2iQz] 



(18) 



They satisfy the boundary conditions T-C{Q,Q) — 1, 
dzT-CiQ, zo) = 0, i/;(0) = 0, and ^WX^o) = 0. The value 
of Zo = (323MeV)~^ is fixed by the rho meson's mass, 
such that Jo(m„zo) = 0, where the lightest mass (labeled 
n — I) corresponds to the rho meson. See also 21, 2^ 

For the soft-wall model ^{z) = k^z'^. The mass 
eigenvalues are given by = 4(n-|-l)K^, with now n = 
corresponding to the rho meson. The integration region 
in Eq. p6|) spans from to infinity. Instead of boundary 
conditions at the upper limit, we require a normalizable 
wave function, J [dz/ z)e~'^ ^ -0^ = 1, and finiteness of H 
at infinity. Using results of H. Grigoryan [j^ . 



H(Q,z) = r(477 + 2)C/(4r;,-l,z2) 
-L^Hn^z% 



k'z' 



1 



(19) 



where lI"' is the generalized Laguerre polynomial and 
U{a,b,'w) is the 2"^^ Kummcr function. 

An analytic form can be obtained for density pQ{b±), 



p+{b^) -JLKo{V2mpb). 

TT 



(20) 



The other densities can be calculated numerically. They 
are shown in Fig. ([2]), along with the charge density [lij . 

Light-front longitudinal densities as well as charge den- 
sities, for hclicity-0 rho-meson are logarithmically di- 
vergent at the origin for both hard-wall and soft-wall 
model, which is evident in ([20]) for the soft-wall model. 
The hard-wall model is more compact than the soft-wall 
model and the helicity-0 rho mesons are more compact 
than the helicity-1 rho mesons. Overall, the distribu- 
tion of longitudinal momentum in position space is more 
compact than that of the charge. 
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FIG. 2: Upper panel: The red solid line is 27rb times P(j^(6), the 
density of helicity-0 p-mesons in the hard- wall AdS/QCD 
model, while the purple dash-dotted line is the corresponding 
result in the soft-wall model. The blue dashed line is 2-nb 
times poib), the charge density of helicity-0 p-mesons in the 
hard-wall model, while the green dash-dot-dotted line is the 
corresponding result in the soft-wall model. Lower panel: the 
same but for pf(b) and pi{b). 



The two independent helicity flip form factors arc 



-l^E. (21) 



However, both B and E vanish in the AdS/QCD model. 

In conclusion, we have studied the distribution within 
extended objects of the matter that carries the compo- 
nent p"*" of the momentum, in a light front viewpoint. 

The examples we used were real nucleons, where we 
used semi-empirical models of the nucleon GPDs as un- 
derlying input, and spin-1 particles, where the underly- 
ing input came from AdS/QCD studies of these states. 
The crucial gravitational form factors can be obtained 
as second moments of the GPDs. There are conceptual 
similarities to the light-front relations of charge distribu- 
tion in the transverse plane to Fourier transforms of the 
electromagnetic form factors. Differences include using 
the gravitational instead of electromagnetic form factors 



and weighting the GPDs with x instead of charge. 

We presented plots that showed the density in the 
entire transverse plane. A qualitative result is that the 
hadrons we study are all more compact when looking 
at the momentum density than when looking at the 
charge (or magnetic) density. We had earlier calculated 
"gravitational radii" from the slope of the gravitational 
form factors obtained for several species of mesons in an 
AdS/QCD model 0, [1^. In addition, we have learned 
that the phenomenon of compactness of the momentum 
distribution and the corresponding smaller root-mean- 
square radius is not limited to mesons which arc studied 
using a purely theoretical AdS/QCD correspondence, but 
is also seen in nucleon distributions based on real data. 

We thank the NSF for support under grant PHY- 
0555600. 



[1] 
[2] 

[3] 
[4] 

[6: 

[7 
[8 
[9 

[lo: 
[11 

[12 
[13 
[14 



[is: 

[16; 
[ir 
[is: 



M. Burkardt, Int. J. Mod. Phys. A 18, 173 (2003) 
[arXiv:hep-ph /0207047 . 

G. A. Miller, P hys. Rev. Lett. 99, 112001 (2007) 
[arXiv:0705.2409l [nucl-th]] . 

C. E. Carlson and M. Vanderhaeghen, Phys. Rev. Lett. 
100, 032004 (2008) |arXiv:0710 .0835 [hep-ph]]. 
X. D. Ji, Phys. Rev. Lett. 78, 610 (1997) 
arXiv:hep-ph/9603249 ; Phy s. Rev. D 55, 7114 
(1997) arXiv:hep-ph /9609381| . 

Z. Abi din and C. E. Ca rlson, Phys. Rev. D 77, 095007 

(2008) garXiv:0801.3839l [hep-ph]]. 

X. D. Ji, J. Phys. G 24, 1181 (1998). 

K. Goeke, M. V. Polyakov and M. Vanderhaeghen, Prog. 

Part. Nucl. Phys. 47, 401 (2001). 

A. V. Belitsky and A. V. Radyushkin, Phys. Rept. 418, 
1 (2005). 

S. Boffi and B. Pa squini, Riv. Nuovo Cim. 30, 387 (2007) 
[arXiv:0711. 26251 [hep-ph]] . 

M. Diehl, Phys. Rept. 388, 41 (2003) 
[arXiv:hep-ph/0307382| . 

M. V. Polyakov, Ph ys. Lett. B 555, 57 (2003) 
farXiv:hep-ph/0 210165| . 

M. Burkardt, Phys. Rev. D 62, 0715 03 (2000) [Erratum- 

ibid. D 66, 119903 (2002)] 'arXiv:hep-ph/0005108 . 

C. E. Carlson and M. Vanderhaeghen, arXiv:080 7T4537l 
[hep-ph]. 

M. Guidal, M. V. Polyakov, A. V. Radyushkin and 
M. Van derhaeghen, Phys. Rev. D 72, 054013 (2005) 
arXiv:hep-ph/0410251 . (For the record, Refs. 0, H 



used slightly different parameters from a preprinted ver- 
sion of this paper, with little difference to the results.) 
L. Frankfurt, C. E. Hyde-Wright, M. Strikman 
and C. Weiss, Phys. Rev. D 75, 054009 (2007) 
^rXiv:hep-ph/0608271 . 



D. Martin, R. G. 
R. S. Thorne, Phys. 



[arXiv:hep-ph/0201127 



Roberts, W. J. 
Lett. B 531, 



Stirling and 
216 (2002) 



A. V. Al'anasev,E7X~Brodsky, C. E. Carlson, Y. C. Chen 
and M. Vanderhaeghen Phys. Rev. D 72, 013008 (2005) 
arXiv:hep-ph/0502013| . 

I. Abidin and C. E. Carlson, Phys. Rev. D 77, 037301 



5 



(2008) [arXiv:0706.26TT] [hep-ph]] . 
[19] Z. Abi din and C. E. Ca rlson, Phys. Rev. D 77, 115021 

(2008) [arXiv:0804.0"2T4l [hep-ph]] . 
[20] J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, Phys. 

Rev. Lett. 95, 261602 (2005) arXiv:hep-ph/0501128 . 
[21] H. R. Grigoryan and A. V. Radyushkin, Phys. Lett. B 

650, 421 (2007) arXiv:hep-ph /0703069; ; Phys. Rev. D 

76, 095007 (20071 |arXiv:0706.1543l [Tiep-ph]]. 



[22] S. J. Brodsky and G. F. de Teramond, Phys. Rev. D 77, 
056007 (2008) arXiv:0707.3859 [hep-ph]].' Phys. Rev. D 
78, 025032 (2008) _arXiv:0804.0452 [hep-ph]]; 

[23] A. Karch, E. Katz, D. T. S on and M. A. St ephanov, Phys. 
Rev. D 74, 015005 (2006) [arXiv:hep-p"h/0602229 . 

[24] H. Grigoryan, unpubhshed. 



